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Abstract. 

Inelastic neutron scattering experiments were performed on single crystals of 
the heavy-fermion compound Celn 3 for temperatures below and above the Neel 
temperature, Tm- In the antifcrromagnetically ordered phase, well-defined spin- 
wave excitations with a bandwidth of 2 meV are observed. The spin waves coexist 
with quasielastic (QE) Kondo-type spin-fluctuations and broadened crystal-field (CF) 
excitations below Tjy. Above TV, only the QE and CF excitations persist, with a weak 
temperature dependence. 
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1. Introduction 

A quantum critical point (QCP) occurs in the heavy fermion (HF) compound Celn 3 
under a hydrostatic pressure of Pc = 25 kbar. This corresponds to a transition at T = 
K from a magnetically ordered to a disordered state. It was found that at least up 
to 23 kbar the ground state is antiferromagnetically ordered with a propagation vector 
k = (1/2,1/2,1/2). At ambient pressure the Neel temperature is TV = 10.1 K and 
the staggered magnetization is m = 0.5 \i& [TJ |2j. Celn 3 becomes superconducting 
below Tc = 200 mK in the vicinity of Pc [HI Ej More generally, the appearance of 
unconventional superconductivity around Pc in HF systems is supposed to be related 
to the enhancement of magnetic fluctuations when TV ~ 0. 

Celn3 crystallizes in the cubic AuCu3 structure, and hence the magnetic 
fluctuations are expected to be three-dimensional (3D). This is confirmed by resistivity 
measurements, which show a Ap = T 3 / 2 non Fermi liquid (NFL) behavior at the critical 
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pressure, as expected for 3D fluctuations in spin- fluctuation theories jSlEl- Recently, 
a new family of HF compounds was discovered [7], Ce m T n In 3m+2n , which is composed 
of alternating layers of Celn3 and TIn 2 stacked along the c-direction with T = Ir, Rh, 
or Co. Their main interest arises from the layered structure that is thought to induce 
2D fluctuations that would enhance the superconductivity more than 3D fluctuations 
near T c [Ej. Indeed, CeCoIn 5 has a superconducting state at ambient pressure below 
Tq = 2.3 K, the highest Tc observed up to now in HF systems 0. These recent results 
motivated us to reinvestigate the common brick of these quasi 2D compounds: the 3D 
heavy fermion Celn 3 . 

Previous inelastic neutron scattering (INS) studies on powder samples of Celn 3 
[TUl ITT] showed the existence of a crystal field (CF) excitation at A C f — 12 meV 
with a large broadening ascribed to the Kondo effect. In this work, we present new 
INS measurements on single crystals, which due to new flux growing methods are now 
sufficiently large for such studies. These measurements show the existence of well-defined 
spin-wave excitations in the magnetically ordered phase, which coexist with quasielastic 
(QE) Kondo-type spin-fluctuations and the previously observed broadened crystal-field 
excitations. 

2. Experimental details 

The single crystals used in this study were grown by the self flux method |12j . A 
first set of neutron scattering measurements was carried out on the time-of-flight 
chopper spectrometer MARI at ISIS, using incident energies of 22 and 60.3 meV. The 
corresponding energy resolution was approximately 2 and 5 meV, respectively. Because 
of the strong neutron absorption cross-section of In, the single crystals were cut to a 
thickness of 2 mm and used in transmission geometry. A total of 30 crystals with total 
weight of 13 g were aligned and glued to a thin aluminium plate, which was mounted 
in a closed-cycle refrigerator. The mosaicity of the assembly was approximately 6°. 
Measurements were performed in the (hhl) plane, with an angle between the incident 
wavevector and the [001] direction of 2.5 and 30°. 

A second set of neutron scattering measurements was made on the triple-axis 
spectrometer IN22 at the ILL (Grenoble, France). Pyrolytic graphite was used for 
the vertically focusing monochromator and the horizontally focusing analyzer. A PG 
filter in the scattered beam was used to suppress higher order contamination, and final 
neutron energies of 8.05 and 14.7 meV were used, giving an energy resolution of 0.5 and 
1.0 meV, respectively. Two thin plates of the assembly used on MARI were aligned on 
an aluminum support and mounted in a helium flow "orange" cryostat. The total mass 
of the sample was about 1 g and the mosaicity 0.5°. The (hhl) scattering plane was 
investigated. 
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3. Results 

3.1. Crystal-field excitation and quasielastic line 

For the two crystals orientations used on MARI and at temperatures of T = 5, 15 
and 150 K, we recorded (Q,E) intensity maps for 0< Q <10 A -1 and E < 30 meV, 
where Q is the wavevector transfer and E is the energy transfer. Despite strong 
phonon contamination from both the Celn3 sample and the Al support, we observed an 
excitation around 12 meV which was independent of the crystal orientation and whose 
intensity decreased with Q. The temperature dependence of the 12 meV excitation 
obtained for Q < 2 A -1 is shown in figure 1. The excitation shows only a very slight 
softening and broadening as the temperature goes from the ordered state at T = 5 K 
to the paramagnetic state at T = 15 K. At much higher temperatures, T = 150 K, 
the excitation is partly suppressed. Figure 2 shows the wavevector dependence of the 
excitation spectrum at T = 5 K. The peak at 12 meV observed at low Q (integration 
over < Q < 2 A -1 ) disappears at higher wave vectors (integration over 7.5 < Q < 7.8 
A -1 ), reflecting its magnetic nature. The contribution around 20 meV, seen only at 
high O/s, is due to phonons. The Q and T dependence of the 12 meV excitation, as 
illustrated in Figs. 1 and 2, clearly show its magnetic origin, and corresponds to the CF 
excitation observed in powder samples by Lawrence et al. ^Oj and Murani et al. jllj . 

The cubic environment of Ce 3+ in Celn3 can be treated in terms of Stevens operators 
by the following CF Hamiltonian: Hqf — -84(04+504) This Hamiltonian splits the 
J = 5/2 spin-orbit level of Ce 3+ into a T 8 quartet and a doublet. In the case of Celn 3 , 
susceptibility measurements indicates that _B 4 > [H], implying that the ground state 
is Tj. This has been confirmed by polarized neutron scattering measurements, which 
also reveals an admixture of r 8 to the ground state T 7 in the ordered phase [T3] . The 
excitation observed here at 12 meV corresponds then to the — ► r§ excitation. 

In addition to the 12 meV CF excitation, it is clear from the data that there 
is an additional quasielastic (QE) component to the scattering. The low-O, data at 
T = 5 K were therefore analyzed using three components: an incoherent elastic peak, 
a quasielastic spin-fluctuation contribution, and a crystal-field excitation. Since the 
different contributions partly overlap, we fixed the width (half width at half maximum, 
HWHM) T QE of the QE contribution to k B T K = 0.86 meV where T K = 10 K is 
the Kondo temperature. The results of the analysis at T = 5 K give a crystal-field 
splitting of Acf = 10.6 ± 0.2 meV and a broadening of the CF levels (HWHM) of 
Tcf — 5.7 ± 0.2 meV. No significant difference in the values of the CF excitations 
were found for the different orientations of the crystal. At T = 15 K, the resulting CF 
splitting is A C f = 9.4 ± 0.2 meV and its width is Tcf — 6.0 ± 0.2 meV. The slight 
softening observed of Acf is related to the splitting of the CF levels below T/y. 
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3.2. Spin wave 

The dispersive mode that appears below Tv was first observed using the MARI chopper 
spectrometer. Because of the 3D character of this mode, the study of its dispersion was 
continued using the IN22 triple-axis spectrometer. Energy scans up to 6 meV energy 
transfer were measured for different values of the wavevector transfer Q. Figure 3 shows 
spectra of Celn 3 at Q = (—1.2,-1.2,0.5) for T = 1.5 and 15 K. This corresponds 
to the reduced wavevector q = Q — r — k = (0.3,0.3,0), where r = (—1,-1,1) 
is a reciprocal lattice vector and k = (—1/2,-1/2,-1/2) is a magnetic propagation 
vector (wavevectors are expressed in reciprocal lattice unit with 1 r.l.u. = 2 it /a, where 
a = 4.689 A is the lattice parameter). Figure 3 clearly shows that the excitation present 
at T < T/v disappears for T > T N . For each wavevector q, the excitation was fitted 
by a Gaussian whose width corresponds to the instrumental resolution. The strong 
dispersion of this excitation in the three main directions of the crystal is shown in figure 
4. Since the direction of the ordered moment in Celn3 is unknown and our sample is in 
a multidomain state, it is not possible to determine the polarization of the excitation; 
it will be assumed to be transverse as for a conventional spin wave. 

Supposing that the CF energy is of the same order as the exchange interaction, 
it seems appropriate to treat the spin wave using the so-called molecular-field random- 
phase approximation ^B|. Writing the exchange Hamiltonian as 

Hex = ^ ] lijJiJj (1) 

gives the general dispersion relation 

E(q) = {[A - M 2 /(q)][A - M 2 /(q+k)]}V2 j (2) 

where k is the propagation vector of the ordered state, /(q) the Fourier transform of 
the exchange interaction 1^-, and A the splitting of the r 7 doublet into r 7j i and r 7)2 
due to the molecular and anisotropy fields. The transition matrix element between the 
ground state r 7) i and the first excited state r 7i2 is given by M 2 = \ < r 7)2 | J a |r 7i i > | 2 
(a = + or -). Considering only the nearest neighbor exchange I\ [along (1,0,0)] and the 
next nearest neighbor exchange J 2 [along (1,1,0)], the exchange interaction is 

/(q) = 2/i [cos (2nq x ) + cos{2irq y ) + cos(2iiq z )} 

+ 2I 2 [cos(2ix(q x - q y )) + cos(2n(q y - q z )) + cos(2n(q z - q x )) 

+ cos(2n(q x + q y )) + cos(2n(q y + q z )) + cos(2n(q z + q x ))}. (3) 

In a two sub-lattice picture constituted of alternating (1,1,1) planes where the spins are 
up or down, Ii would correspond to the antiferromagnetic exchange between spins from 
different sublattices and I2 would correspond to the ferromagnetic exchange between 
spins from the same sublattice. A global fit to the dispersion curves of equation 
(J2J) is shown in figure 4. It gives the exchange terms M 2 I\ = —0.354 ± 0.040 meV 
and M 2 I 2 = 0.028 ± 0.013 meV and the T 7 splitting A = 2.812 ± 0.480 meV. The 
corresponding spin- wave gap is : 5 = E(0) = 1.28 meV. 
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The exchange energy can be expressed by E ex = 2J ^ +1) 7(0). Taking the wave 
functions of r 71 and T 72 when there is neither exchange nor distortion ^7] gives us an 
approximated square transition matrix element M 2 = 2.78 thus 1(0) = —0.64 ± 0.18 
meV. With J = 5/2, we obtain E ex ~ 3.71 meV. The ratio A CE /\E ex \ ~ 3 permits 
then to conclude that the CF energy is of the same order than the exchange interaction 
that justifies the treatment of the spin wave using the molecular-field random-phase 
approximation. 

4. Discussion 

4-1. Kondo temperature 

We have analyzed the Q-independent excitations using two contributions, one 
corresponding to quasielastic scattering and the other to a crystal-field transition. Since 
these excitations overlap with each other and with the incoherent peak (see figure 2), the 
energy, width, and weight of the different contributions extracted from the data depend 
on the model used for the fitting. This is illustrated by the scatter in the estimates 
of A CF , Tcf, and Tq E found in the literature [TUJ HU HH1 UH] • At low temperatures, 
the QE excitation corresponds to the spin fluctuations of the T 7 electronic level, which 
are attributed to the Kondo effect. Its width (HWHM) Tq E corresponds to the Kondo 
temperature. A Kondo temperature of Tk = 10 K at a pressure of P = 19 kbar was 
obtained in 115 In-NQR measurements by Kawasaki et al. J20 3 . It is likely that below this 
pressure the anomaly at TV in the 1/Tl temperature variation masks the Tk anomaly. 
We assume therefore that Tk has a very slight variation for pressures below P = 19 
kbar, so that it is only slightly smaller than 10 K at ambient pressure. Resistivity 
measurements give a Kondo scattering maximum at Tm ~ 50 K |3], which corresponds 
to the high Kondo temperature T K = ^T K (A CF /k B ) 2 [HI INS gives A CF = 10.6 
meV which with T K = 50 K leads to Tk — 8 K. It is then reasonable to approximate 
the Kondo temperature by Tk — 10 K in Celn3 at ambient pressure. For the INS 
measurement carried out at T = 5 K < Tk, we choose consequently to fix the QE 
width to Tqe = ^bTk = 0.86meV. The fact that Tcf = 5.7 meV is relatively large 
compared to Tq E = 0.86 meV is due to the larger spin fluctuations associated with the 
Tg first excited level. In fact, the CF transition T 7 Tg has a width corresponding to 
the convolution of the T 7 and r 8 widths T]j E = Tq E and T^ E . We can approximate 
T C f — ^ T qe + Fqe and thus estimate that ~ T C f — ^qe = 4.84 meV, which is 
effectively 5 times larger than Tq E . 

4-2. Anisotropy 

Even if the CF level is broad, it seems that the propagating spin waves can be described 
by considering T 7 as an "isolated" doublet split by the anisotropy and exchange fields. 
Since this doublet originates from a sixfold degenerate J = 5/2 level, an anisotropy 
gap is expected in the spin-wave spectrum. We recall that in the isotropic case, the 
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spin-wave dispersion has no gap, since only the exchange is responsible for the splitting 
A. 

The gap 5 of the spin wave excitation corresponds to its minimum energy that 
is obtained at q = (0,0,0). Since this excitation overlaps with the incoherent 
elastic scattering that reflects the experimental resolution and with the quasielastic 
contribution of the spin fluctuations, the gap was not correctly obtained in the present 
neutron experiment. The given value of S corresponds to the best fit to the data using the 
spin-wave model described above. An attempt to directly observe the gap using a smaller 
incident energy (and hence better resolution) on the cold triple-axis spectrometer IN12 
(ILL) failed because of the strong indium absorption cross-section. A determination of 
the energy gap was also made by low-temperature specific-heat measurements which 
exhibited a spin wave contribution that leads to a gap 5 = 0.69 meV [22]. This suggests 
that the gap energy of 1.28 meV deduced from our INS measurements is somewhat 
overestimated. 

The result of our spin-wave analysis is in good agreement with calculations made 
by Wang et al. for Ce 3+ ions under cubic crystal field and effective exchange field 
|17j . They consider a molecular field along the c-axis, which can be approximated 
by g\i B H m = k B T N , and calculate the level splitting as a function of g[i B H m j B^. 
Knowing that for Celn 3 T N = 10.1 K, A CF = 10.6 meV and hence B 4 = 2.94 10" 2 
meV (A C f = 360 _B 4 ), we have g^ B H m = 8.71 10 _1 meV ~ 30 B A . For this value of 
gfi B H m / B^, Wang et al. calculate a T 7 splitting of about 40£> 4 ~ 1.2 meV. Although 
we suppose that in our case the molecular field is not along the c-axis (see below), we 
believe that Wang's calculations give the right order of the splitting A. 

There are several indications for that the moments are not along the c-axis. If 
the CF alone imposes the easy magnetization axis, the positive sign of B^ in the CF 
Hamiltonian implies that the moments are along the [1,1,1] direction in the ordered 
phase. This is also supported by nuclear quadrupolar resonance measurements |23j. 
Although it is difficult to determine the moment direction of a cubic compound by 
neutron diffraction, some recent measurements also suggest that the moments are not 
along the c-axis [2*2"] . 

4-3. Comparison with CePd 2 Si2 

Celn3 seems to have a behavior very similar to the HF compound CePd2Si2 which is 
antiferromagnetically ordered at ambient pressure and has the same kind of excitation 
spectrum with coexisting QE, CF, and spin-wave excitations [21]. The characteristic 
energies are very close and, moreover, unconventional superconductivity appears at a 
similar external pressure [3]. A summary of the characteristic physical quantities of 
those two compounds is given in Table 1. A focus on the few differences between those 
two compounds can be made as following. Firstly CePd2Si2 is tetragonal while Celn 3 is 
cubic. Some characteristics such as the CF splitting or the susceptibility anisotropy can 
be directly linked to anisotropy effects that are resulting from those lattice structures. 
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Table 1. Characteristic physical quantities of Celri3 and CePd2Si2. Data are taken 

from 0, g], El, Ha, ESI, H3, ESI- 





Celn 3 


CePd 2 Si 2 




10.1 K 


10 K 


k 


(1/2,1/2,1/2) 


(1/2,1/2,0) 


m 


0.50 fi B 


0.62 fi B 




0.69 meV 


0.83 meV 


Bandwidth 


~ 2 meV 


~ 2 meV 


0p 


-50 K 


-47 K * 


T K 


10 K 


10 K 


A CF 


123 K 


220 and 280 K 


Pc 


26.5 kbar 


28.6 kbar 


T c 


200 mK 


430 mK 




130 mJ.moH.K" 2 


250 mJ.moH.K- 2 



* : DC susceptibility measurement along a-axis 
** : 7 = C/T at low T 

Then, the characteristic quantities Tc and 7 are twice bigger in CePd2Si2 than in Celn 3 , 
that could traduce an enhancement of the spin fluctuations in the tetragonal CePd2Si2- 
As a supplement to the studies of quasi 2D Ce m T n In 3m+ 2 n compounds where Tc can 
reach more than 2 K, an alternative way to study the influence of anisotropy in strongly 
correlated systems could consequently remain in a detailed comparison between cubic 
Celn 3 and tetragonal CePd 2 Si 2 . 

5. Conclusion 

We have performed the first inelastic neutron scattering measurements on single crystals 
of Celn 3 . The combined use of time-of-flight and triple- axis spectrometers allowed us to 
study both local and dispersive low-energy excitations. We found the coexistence below 
Tjy of a spin-wave with a quasielastic component and a crystal-field excitation. The CF 
levels are substantially broadened by Kondo-type spin fluctuations, which indicates the 
proximity to an intermediate valence state where the CF collapse. Despite the smearing 
of the CF levels, the observation of spin waves means that the ground state level in Celn 3 
is also relatively well-defined. The challenge is then to understand how these excitations 
evolve on going toward a QCP where the hierarchy between the Kondo effect and the 
RKKY coupling is reversed while the CF should be close to a collapse. 
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Figure 1. Energy scans at different temperatures obtained on MARI spectrometer 
using an incident energy of 60.3 meV and integrating over Q values between and 2 
A -1 . The solid line shows the total fit for T = 5 K. 
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Figure 2. Energy scans at different wavevectors obtained on the MARI spectrometer 
using an incident energy of 60.3 meV at T — 5 K. The long dashed, short dashed, and 
dotted lines show the CF, QE, and incoherent contributions of the small Q spectrum, 
respectively, while the solid line shows the total fit. 
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Figure 3. Energy scans at Q = (—1.2, —1.2, 0.5) for temperatures of T = 1.5 and 15 
K obtained on the IN22 spectrometer using a final energy of 14.7 meV. 
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Figure 4. Spin-wave dispersion of Celn3 along the three main directions (0,0, q), 
(q,q,0) and (q,q,q) at T = 1.5 K, obtained on the IN22 spectrometer using a final 
energy of 14.7 mcV. The solid line is a fit obtained using the molecular-field random- 
phase approximation. 



